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Abstract. We present model spectra of stellar populations with variable chemical composition. 
We derived the [a/Fe] abundance ratio of the stars of the most important libraries (ELODIE, 
CFLIB and MILES) using full spectrum fitting and we generated PEGASE.HR models resolved 
in [a/Fe]. We used a semi-empirical approach that combines the observed spectra with synthetic 
stellar spectra. We tested the models using them to derive [a/Fe] in galaxies and star clusters 
using full spectrum fitting. The present models are available from http : //ulyss . univ-lyonl . f r 
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1. Introduction 
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^""^ , The comparison between models and observations is the key to derive the metallicity, 
■ age and star formation history of galaxies and star clusters. The evolutionary synthesis of 
galaxies consists in choosing an initial mass function (IMF) and following the evolution 
of the population over the time according to some theoretical evolutionary tracks. This 
+3 [ approach predicts the distribution of the stars in the effective temperature (T^ff) vs. 

surface gravity (log g) plan at any epoch, and the final step is to draw spectra from a 
i— —i | library to assemble the integrated spectrum. 

The m odelling of stel lar populations steadily progressed since the first attempts in the 
80s (e.g. B ruzuall 119831) . Beside the metallicity, age and history, the detailed chemical 
• composition became an important parameter to measure. This allows one to trace the 
y—{ \ source of metal enrichment, because of the different yields of type la and II supernovae. 

In turn, since both SN types correspond to different time-scales, this informs about the 
duration of the enrichment process. A population that produced its metal and formed 
£N) ■ its stars before the SNia onset will be deficient in Fe (or as more commonly said, over- 
abundant in a-elements, like O and Mg). To vary the yield in the models, both the 
stellar library and the evolutionary tracks have in principle to be changed. However, 



Potter et al.l (|2007l) have shown that the tracks essentially depend on the total metallicity, 
rather than of the actual mixture. Therefore, the scaled-solar tracks can reasonably be 
used for any composition, provided they are associated with stellar spectra of the same 
j_j • total metallicity. 

The models predict colours, used for SED fitting, spectra, and line-strength indices. 
Indices are equivalent-widths defined for a number of prominent spectral features and 
their study brings an approximate handle on the physical parameters (age, metallicity 
and/or detailed composition). 

Several authors modeled indices with variable abundances of a-elements relative to Fe 



(|Weiss et al.lll995t iTrager et aljliooot iThomas et al.l [2003). The approach consisted in 



evaluati ng the response of ind ividual indices to changes of composition using theoretical 
spectra (|Tripicco fe Belllll995l) . Then, assuming that the abundance pattern in the library 
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is matching that of the solar neighborhood (justified by the fact that the observed libraries 
consist of nearby stars), the indices could be estimated for any T e g log g, [Fe/H] and 
[a/Fe]. In the visible region of the spectrum, the variation of [a/Fe] results in a prominent 
change of the Mg^ index near 5175 A, and therefore a combination of this index with 
some Fe indices could be used to derive the [a/Fe] in galaxies and clusters. 

Models of indices were developing faster and further than the spectral predictions, 
because of the lack of stellar libraries calibrated in relative flux. However, the quality 
of the libraries also progressed. Libraries of about 1000 stars covering a wide range 
of the parameter space, with an extended coverage of the optical spectrum at mediun 
spectral resolu tion and a fair flux cali bration became available. The three most important 
are ELODIE (IPrugniel et alJ l2007bL R=A/AA «10000, AA = 3900 - 6800 A), CFLIB 



( |Valdes et al 
120061 



2004 Rss4000, A A = 3460 
2000, AA = 3536- 7410 A). 



9400 A) and MILES (jSanchez-Blazquez et al 



In this paper, we present the first spectral models with variable [a/Fe] based on these 
libraries and on precise measurements of [a/Fe] in their stars. 



2. Method 

Previ ous attempts to produce semi -empirical spectral mo dels with variable [a/Fe] were 



made in IPrugniel et all (|2007al) . and IWalcher et all (|2009h . The latter one uses popula- 



tion models based on synthetic spectra to predict the differential change of each spectral 
element when [a/Fe] is varied. This differential effect is then used to 'correct' empirical 
models to the desired a-elements enrichment. The former method determines the differ- 
ential changes of the stellar spectra and 'correct' the library before computing models at 
different [a/Fe]. Both approachs are equivalent, and are based on the assumption that 
the libraries follow the abundance pattern of the solar neighborhood (i.e. they adopt 
a functional dependence of [a/Fe] as a function of [Fe/H]). Models of Lick indices also 
make the same assumption. 

Although the assumption of a relation between [a/Fe] and [Fe/H ] is statistically ver - 
ified, [a/Fe] present a significant dispersion at a given [Fe/H] (e.g. iMilone et al.ll201lf ). 
When the library is interpolated to extract the spectra entering in the synthesis, the 
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Figure 1. Distribution of [a/Fe] in the three l ibrarie s, compared to the distribution in 

IMilone et aP lj2011f ). 



Population models with a/Fe 
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stochastic effects can become significant. For example, for a local interpolation only a 
small number of library stars are used near a given location in the T e g, log g and [Fe/H] 
space, and their mean [a/Fe] may be biased. In order to avoid this effect, the improve- 
ment introduced in the present work is to determine [a/Fe] in the individual stars of 
the library. Each stellar spectrum is then corrected to a set of given [a/Fe] applying 
differential effects inferred from a theoretical library. An alternative may have been to 
measure [a/Fe] and directly interpolate the spectra in the four-dimensional space. This 
would have avoided the dependency toward a theoretical library, but required some ex- 
trapolations, because of the non-uniform distribution of the stars in the [Fe/H] vs. [a/Fe] 
plan,, 

The approach is finally the following: 

- Measure [a/Fe] in the spectra of the libraries. 

- Correct the library to some given values of [a/Fe] using differential effects computed 
in a theoretical library. 

- Build interpolators for these corrected libraries, (these are functions returning a 
spectrum for given T e g, log g, and [Fe/H]. 

- Make model predictions at the chosen [a/Fe]. 



3. Determination of [a/Fe] in the stars 

We used full-spectrum fitting with the ULySS package ( Koleva et al. 20091 ) to measure 
[a/Fe]. Full-spectrum fit ting was successful to measure the other atmospheric parameter s 
in the various libraries (|Wu et al.ll201ll IPrugniel et al. 2011 : Koleva fc Vazdekisl 2011 ). 
and the precision and reliability was found consistent with detailed analysis using high- 
resolution spectroscopy. 

We measured [a/Fe] by comparing the ELODIE spectra to the grid of ICoelho et al 



( 2005h . We minimized the X 1 residuals to derive the four atmospheric parameters 



We 



measured [a/Fe] only for 3500 < T e // < 7000, because of the limits of the grid of 
synthetic spectra. We assumed that the remaining stars follow the solar neighborhood 
abundance pattern. Then, we used the ELODIE library as a reference to measure [a/Fe] 
in the other libraries, using the same full-spectrum fitting method. 

The measurements c ompa re well with the most precise measurements from the compi- 
lation of iMilone et al.l (|201lf ). The distr ibution of [a/Fe] vs . [Fe/H] is presented in Fig.Q] 



for the three libraries and compared to IMilone et al. (2011) 



4. Population models resolved in [a/Fe] 

We generated pairs of interpopators for [a/Fe] £ {0,0.4} d ex, and computed grid s 
of single stellar populations (SSPs) models with PEGASE.HR (|Le Borgne et al.l 12004 ). 
These pairs of grids were then grouped together to make four dimensional arrays (wave- 
length, age, [Fe/H] and [a/Fe]) that can be used with ULySS to determine the three 
SSP-equivalent parameters. Figure [5] is a snapshot of models made with the three li- 
braries. 



We used these models to analyze Galactic globular clusters spectra as in iKoleva et al 



( 20081) . and we found results consistent with the determinations based on sp ectroscopy of 



resolv ed stars. We also applied the method to study spectra of galaxies from IKoleva et al 



( 20111 ) and we found that the results are consistent with those from Lick indices, but they 
are more precise. The gain in precision reflects the better usage of the information. Mgf, 
accounts for less than 20 % of the total opacity change in the optical range (A > 4500 A) 
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when [a/Fe] varies, and using the whole spectral range with full-spectrum fitting is a 
clear advantage. 



5. Prospective 

This semi-empirical method, used to build models of stellar populations resolved in 
[a/Fe], can straightforwardly be applied to resolve other individual elements. The key 
poin t is to have r e liable synthetic grids of spectra to determine the differential effects. 
The lCoelho et al.l (|2005f ) grid already resolve Ca independently of the a-elements, and 



other grids, varying other elements will be computed in the future. 



The models presented here are available at http://ulyss.uly-lyonl.fr They can 



be used with the ULySS package to measure [a/Fe] in galaxies and clusters. 



References 

Bruzual, G. 1983, ApJ, 273, 105 

Coelho, P., Barbuy, B., Melendez, J., Schiavon, R. P., & Castilho, B. V. 2005, A&A, 443, 735 
Dotter, A., Chaboyer, B., Ferguson, J. W., et al. 2007, ApJ, 666, 403 

Koleva, M., Gupta, R., Prugniel, P., & Singh, H. 2008, in Astronomical Society of the Pacific 
Conference Series, Vol. 390, Pathways Through an Eclectic Universe, ed. J. H. Knapen, 
T. J. Mahoney, & A. Vazdekis, 302 
Koleva, M., Prugniel, P., Bouchard, A., & Wu, Y. 2009, A&A, 501, 1269 
Koleva, M., Prugniel, P., de Rijcke, S., & Zeilinger, W. W. 2011, MNRAS, 417, 1643 
Koleva, M. & Vazdekis, A. 2011, ArXiv e-prints, 1111.5449 
Le Borgne, D., Rocca-Volmerange, B., Prugniel, P., et al. 2004, A&A, 425, 881 
Milone, A. D. C, Sansom, A. E., & Sanchez-Blazquez, P. 2011, MNRAS, 414, 1227 
Prugniel, P., Koleva, M., Ocvirk, P., Le Borgne, D., & Soubiran, C. 2007a, in IAU Symposium, 

Vol. 241, IAU Symposium, ed. A. Vazdekis & R. F. Peletier, 68-72 
Prugniel, P., Soubiran, C, Koleva, M., & Le Borgne, D. 2007b, arXiv:astro-ph/0703658 
Prugniel, P., Vaughn, I., & Koleva, M. 2011, A&A, 531, A165 

Sanchez-Blazquez, P., Peletier, R. F., Jimenez- Vicente, J., et al. 2006, MNRAS, 371, 703 
Thomas, D., Maraston, C, & Bender, R. 2003, MNRAS, 339, 897 
Trager, S. C, Faber, S. M., Worthey, G., & Gonzalez, J. J. 2000, A3, 119, 1645 
Tripicco, M. J. & Bell, R. A. 1995, AJ, 110, 3035 




4000 5000 6000 7000 8000 9000 

Wavelength, A 

Figure 2. Single stellar population models with PEGASE.HR and ELODIE, MILES and 
CFLIB. The age is eight Gyr and the metallicity is solar. 
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